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A purely algebraic algorithm for computation of invariants (generalized Casimir 
operators) of Lie algebras by means of moving frames is discussed. Results on 
i-S^ . the application of the method to computation of invariants of low-dimensional 

Lie algebras and series of solvable Lie algebras restricted only by a required 
(~| . structure of the nilradical are reviewed. 

a 

1 Introduction 

! The invariants of Lie algebras are one of their defining characteristics. They have 

\^ I numerous applications in different fields of mathematics and physics, in which 

■ Lie algebras arise (representation theory, integrability of Hamiltonian differential 
, equations, quantum numbers etc). In particular, the polynomial invariants of 
I a Lie algebra exhaust its set of Casimir operators, i.e., the center of its universal 
' enveloping algebra. This is why non-polynomial invariants are also called gener- 
. alized Casimir operators, and the usual Casimir operators are seen as 'specific' 

generalized Casimir operators. Since the structure of invariants strongly depends 

■ on the structure of the algebra and the classification of all (finite-dimensional) 
. Lie algebras is an inherently difficult problem (actually unsolvabl^), it seems to 

5^ I be impossible to elaborate a complete theory for generalized Casimir operators 

in the general case. Moreover, if the classification of a class of Lie algebras is 
known, then the invariants of such algebras can be described exhaustively. These 
problems have already been solved for the semi-simple and low-dimensional Lie 
algebras, and also for the physically relevant Lie algebras of fixed dimensions. 

The standard method of construction of generalized Casimir operators consists 
of integration of overdetermined systems of first-order linear partial differential 



^The problem of classification of Lie algebras is wild since it includes, as a subproblem, the 
problem on reduction of pairs of matrices to a canonical form [10]. For a detailed review on 
classification of Lie algebras we refer to [17]. 
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equations. It turns out to be rather cumbersome calculations, once the dimension 
of Lie algebra is not one of the lowest few. Alternative methods use matrix rep- 
resentations of Lie algebras. They are not much easier and are valid for a limited 
class of representations. 

In our recent papers [3-7] we have developed the purely algebraic algorithm for 
computation of invariants (generalized Casimir operators) of Lie algebras. The 
suggested approach is simpler and generally valid. It extends to our problem the 
exploitation of the Cartan's method of moving frames in Fels-Olver version [9]. 
(For modern development of the moving frames method and more references see 
also [14,15].) 

2 Preliminaries 

Consider a Lie algebra g of dimension dimg = n < oo over the complex or real 
field F (either F = C or F = M) and the corresponding connected Lie group G. 
Let 0* be the dual space of the vector space g. The map Ad* : G GL(g*) defined 
for any (7 € G by the relation 

(Ad*a;,u) = {x, Adg-iu) for all x E g* and n € g 

is called the coadjoint representation of the Lie group G. Here Ad: G — > GL(g) 
is the usual adjoint representation of G in g, and the image Ada of G under Ad 
is the inner automorphism group Int(g) of the Lie algebra g. The image of G 
under Ad* is a subgroup of GL(g*) and is denoted by Ad^. 

The maximal dimension of orbits of Adg is called the rank of the coadjoint 
representation of G (and g) and denoted by rankAd^. It is a basis independent 
characteristic of the algebra g. Orbits of this dimension are called regular ones. 

A function F G C°°(r2), where Q is a domain in g*, is called a (global in fl) 
invariant of Adg if F(Ad*a;) = F{x) for all g G G and x G O, such that Ad*a; € O. 
The set of invariants of Adg on Q is denoted by Inv(Ad^) without an explicit 
indication of the domain Q. Let below O is a neighborhood of a point from 
a regular orbit. It can always be chosen in such a way that the group Ad^ 
acts regularly on Q. Then the maximal number A^^ of functionally independent 
invariants in Inv(Ad^) coincides with the codimension of the regular orbits of Ad^, 
i.e., it is given by the difference Ng = dimg — rankAd^. 

To calculate the invariants explicitly, one should fix a basis £ = {ei, . . . , e„} 
of the algebra g. It leads to fixing the dual basis £* = {e^, . . . ,e*} in the dual 
space g* and to the identification of Int(g) and Ad^ with the associated matrix 
groups. The basis elements ei, . . . , e„ satisfy the commutation relations [e^, ej] = 
c^jCk, where c^^- are components of the tensor of structure constants of g in the 
basis £. Here and in what follows the indices i, j and k run from 1 to n and the 
summation convention over repeated indices is used. Let x — > x = (xi, . . . , x„) be 
the coordinates in g* associated with £*. 

It is well known that there exists a bijcction between elements of the universal 
enveloping algebra (i.e., Casimir operators) of g and polynomial invariants of g 
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(which can be assumed defined globally on q*). See, e.g., [1]. Such a bijection is 
established, e.g., by the symmetrization operator Sym which acts on monomials 
by the formula 

Sym(eii • • • J = ^ ^ ei^, • • • Cj,, , 

' aeSr 

where take values from 1 to n, r G N. The symbol Sr denotes the 

permutation group consisting of r elements. The symmetrization also can be 
correctly defined for rational invariants [1]. If Int(Ad^) has no a functional basis 
consisting of only rational invariants, the correctness of the symmetrization needs 
an additional investigation for each fixed algebra g since general results on this 
subject do not exist. After symmetrized, elements from Int(Adgr) are naturally 
called invariants or generalized Casimir operators of g. The set of invariants of g 
is denoted by Inv(g). 

Functionally independent invariants F^xi, . . . ,Xn), I = I,. . . ,Ng, forms a 
functional basis {fundamental invariant) of Inv(AdQ) since any element from 
Inv(Ad^) can be (uniquely) represented as a function of these invariants. Ac- 
cordingly the set of SymF'(ei, . . . , e„), 1 = 1,..., Ng, is called a basis of Inv(g). 

In framework of the infinitesimal approach any invariant F{xi, . . . , Xn) of Ad^ 
is a solution of the linear system of first-order partial differential equations [1,2,16] 
XiF = 0, i.e., c^jXkFxj = 0, where Xi = cfjXkdxj is the infinitesimal generator 
of the one-parameter group {AdQ{exp eei)} corresponding to e^. The mapping 
Cj Xi gives a representation of the Lie algebra g. 

3 The algorithm 

Let Q = AdQ X g* denote the trivial left principal Ad^-bundle over g*. The right 
regularization R of the coadjoint action of G on g* is the diagonal action of Adg 
on ^ = Adg. X g*. It is provided by the map 

Rg{Adl, x) = {Adl ■ Ad*_i , Ad*gx), g,heG, xeq*. 

The action R on the bundle Q = Ad^ x g* is regular and free. We call R the 

lifted coadjoint action of G. It projects back to the coadjoint action on g* via 
the Adg-equivariant projection ttq*: G 9*- Any lifted invariant of Adg. is 
a (locally defined) smooth function from ^ to a manifold, which is invariant with 
respect to the lifted coadjoint action of G. The function I: Q ^ q* given by 
I = I{Ad*,x) = Ad*a; is the fundamental lifted invariant of Ad^, i.e., I is 
a lifted invariant and any lifted invariant can be locally written as a function of I 
in a unique way. Using an arbitrary function F{x) on g*, we can produce the 
lifted invariant F o I of Ad^ by replacing x with I = Ad*x in the expression 
for F. Ordinary invariants arc particular cases of lifted invariants, where one 
identifies any invariant formed as its composition with the standard projection vTg*. 
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Therefore, ordinary invariants are particular functional combinations of lifted ones 
that happen to be independent of the group parameters of Ad-Q. 

The essence of the normalization procedure by Fels and Olver can be presented 
in the form of on the following statement. 

Proposition 1. Suppose that X = (Ti, . . . is a fundamental lifted invariant, 
for the lifted invariants Zj^, . . . , Zj^ and some constants ci, . . . , Cp the system 
Zj^ = ci , . . . , = Cp is solvable with respect to the parameters Ok^ , . . . , 9^^ 
and substitution of the found values of Ok^, ■ ■ ■ , Okp into the other lifted invariants 
results in m = n — p expressions Zi, I = 1,. . . ,m, depending only on x's. Then 
p = rank Ad^, m = Ng and Zi, . . . , Zm form a basis of Inv(Ad^). 

The algebraic algorithm for finding invariants of the Lie algebra g is briefly 
formulated in the following four steps. 

1. Construction of the generic matrix B{6) of Ad^. B{9) is the matrix of 
an inner automorphism of the Lie algebra g in the given basis ei, . . . , e^, 9 = 
{9i, . . . ,9r) is a complete tuple of group parameters (coordinates) of Int(g), and 
r = dimAdf; = dimlnt(g) = n — dimZ(g), where Z(g) is the center of g. 

2. Representation of the fundamental lifted invariant. The explicit form of 
the fundamental lifted invariant Z = (Zi, . . . ,Zn) of Ad^ in the chosen coordi- 
nates {9, x) in Ad^ x g* is X = x ■ B{9), i.e., 

(Ti, . . . ,Zn) = (xi, ...,x„)- B{9i, . . .,9r). 

3. Elimination of parameters by normalization. We choose the maximum pos- 
sible number p of lifted invariants Xj^, . . . , Zj^, constants ci, . . . , Cp and group 
parameters 9^1, ... , 9^^ such that the equations Zj-^ = ci , . . . , Zj^ = Cp are solvable 
with respect to 9k^, . . . , 9kp. After substituting the found values of 0^^, . . . , 9kp 
into the other lifted invariants, we obtain Ng = n — p expressions F\xi, . . . , Xn) 
without 0's. 

4. Symmetrization. The functions F'^{xl, . . . ,Xn) necessarily form a basis of 
Inv(AdQ). They are symmetrized to SymF'(ei, . . . , e„). It is the desired basis 
of Inv(g). 

Our experience on the calculation of invariants of a wide range of Lie algebras 
shows that the version of the algebraic method, which is based on Proposition [H 
is most effective. In particular, it provides finding the cardinality of the invariant 
basis in the process of construction of the invariants. The algorithm can in fact 
involve different kinds of coordinate in the inner automorphism groups (the first 
canonical, the second canonical or special one) and different techniques of elim- 
ination of parameters (empiric techniques, with additional combining of lifted 
invariants, using a floating system of normalization equations etc). 

Let us underline that the search of invariants of a Lie algebra g, which has been 
done by solving a linear system of first-order partial differential equations under 
the conventional infinitesimal approach, is replaced here by the construction of 
the matrix B(6) of inner automorphisms and by excluding the parameters 9 from 
the fundamental lifted invariant Z = x ■ B[9) in some way. 
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4 Illustrative example 

The four-dimensional solvable Lie algebra g has the following non-zero commu- 
tation relations 

[62,63] =ei, [61,64] = (1 + 6)61, [62,64] =62, [63,64] = 663, |6| < 1. 

Its nilradical is three-dimensional and isomorphic to the Weil-Heisenberg alge- 
bra 93.1. (Here we use the notations of low-dimensional Lie algebras according to 
Mubarakzyanov's classification [11].) 

We construct a presentation of the inner automorphism matrix B{9) of the Lie 
algebra 5, involving second canonical coordinates on Adc as group parameters 6. 
The matrices ade^, i = 1, . . . , 4, of the adjoint representation of the basis elements 
ei, 62, 63 and 64 respectively have the form 
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The inner automorphisms of 34 § are then described by the triangular matrix 

3 

B{9) = exp(6'iadej • exp(-6'4ade4) 

i=l 

/g(i+fe)e4 _^3ee4 Q^^be^ ^,^2^3 + (1 + 5)^ a 
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Therefore, a functional basis of lifted invariants is formed by 
Ti = e(i+^)^4xi, 

l3 = e'''^{e2Xi+X3), 

T4 = {be293 + (1 + b)ei)xi + 6*2X2 + 66*3X3 + X4. 

Further the cases 6 = — 1 and b ^ —1 should be considered separately. 

There are no invariants in case b ^ —1 since in view of Proposition[T]the number 
of functionally independent invariants is equal to zero. Indeed, the system Zi = 1, 
X2 = Z3 = 2^4 = is solvable with respect to the whole set of the parameters 9. 
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It is obvious that in the case 6 = — 1 the element ei generating the center Z(gjg) 
is an invariant. (The corresponding Ufted invariant T\ = x\ does not depend on 
the parameters 9.) Another invariant is easily found via combining the lifted 
invariants: X1X4 — 12^2, = 3^13^4 — 3^23^3- After the symmetrization procedure we 
obtain the following polynomial basis of the invariant set of this algebra 

6263 + 6362 
ei, 6164 2 ■ 

The second basis invariant can be also constructed by the normalization technique. 
We solve the equations Z2 = 2^3 = with respect to the parameters 62 and ^3 and 
substitute the expressions for them into the lifted invariant 24. The obtained 
expression X4 — x^xzjxx does not contain the parameters Q and, therefore, is 
an invariant of the coadjoint representation. For the basis of invariants to be 
polynomial, we multiply this invariant by the invariant x\. It is the technique 
that is applied below for the general case of the Lie algebras under consideration. 

Note that in the above example the symmetrization procedure can be assumed 
trivial since the symmetrized invariant 6164 — 5(6263 + 6362) differs from the non- 
symmetrized version 6164 — 6263 (resp. 6164 — 6362) on the invariant ^ei (resp. 
— \e\). If we take the rational invariant 64 — 6263/61 (resp. 64 — 6362/61), the 
symmetrization is equivalent to the addition of the constant \ (resp. — 5)- 

Invariants of 54 § were first described in [16] within the framework of the in- 
finitesimal approach. 

5 Review of obtained results 

Using the moving frames approach, we recalculated invariant bases and, in a 
number of cases, enhanced their representation for the following Lie algebras (in 
additional brackets we cite the papers where invariants bases of the same algebras 
were computed by the infinitesimal method): 

• the complex and real Lie algebras up to dimension 6 [3] ( [8, 12, 16]); 

• the complex and real Lie algebras with Abelian nilradicals of codimension 
one [4] ( [18]); 

• the complex indecomposable solvable Lie algebras with the nilradicals iso- 
morphic to ^5q , n = 3, 4, . . . (the nonzero commutation relations between 
the basis elements ei, e„ of Zq are exhausted by [6fc,6„] = 6fe-i, 
A: = 2,...,n-1) [4] ( [13]); 

• the nilpotent Lie algebra io{n) oi n x n strictly upper triangular matri- 
ces [4,5] ( [20]); 

• the solvable Lie algebra t(n) of n x upper triangular matrices and the 
solvable Lie algebras st(n) of n x n special upper triangular matrices [5-7] 
( [20]); 

• the solvable Lie algebras with nilradicals isomorphic to iQ{n) and diagonal 
nilindependent elements [5-7] ( [20]). 
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Note that earlier only conjectures on invariants of two latter families of Lie 
algebras were known. Moreover, for the last family the conjecture was formulated 
only for the particular case of a single nilindependent element. Here we present 
the exhaustive statement on invariants of this series of Lie algebras, obtained 
in [7]. 

Consider the solvable Lie algebra ty{n) with the nilradical NR(ty(n)) isomor- 
phic to to(n) and s nilindependent element fp,p= 1, . . . , s, which act on elements 
of the nilradical in the way as the diagonal matrices Tp = diag(7pi, . . . , 7p„) act on 
strictly triangular matrices. The matrices Tp, p = 1, . . . , s, and the unity matrix 
are linear independent since otherwise NR(tY(n)) ^ to(n). The parameter matrix 
7 — ilpi) defined up to nonsingular s x s matrix multiplier and homogeneous 
shift in rows. In other words, the algebras ty(n) and ty'(n) are isomorphic iff there 
exist A G Ms,s(F), det A ^ 0, and /x G such that 

s 

I'pi = XI '^pp'^P'i + l^p^ p=l,...,s, i = l,...,n. 
p'=i 

The parameter matrix 7 and 7' are assumed equivalent. Up to the equivalence 

the additional condition Tr Fp = = can be imposed on the algebra 

parameters. Therefore, the algebra ty(n) is naturally embedded into st(n) as a 
(mega)ideal under identification of NR(ty(n)) with to(n) and of fp with Tp. 

We choose the union of the canonical basis of NR(ty(n)) and the s-clement 
set {fp,p = 1, . . . , s} as the canonical basis of ty(n). In the basis of NR(ty(n)) we 
use 'matrix' enumeration of basis elements Cij, i < j, with the 'increasing' pair 
of indices similarly to the canonical basis {E^j, i < j} of the isomorphic matrix 
algebra to(n). 

Hereafter E^j (for the fixed values i and j) denotes the n x n matrix {Sn'Sjj') 
with i' and j' running the numbers of rows and column correspondingly, i.e., the 
nxn matrix with the unit on the cross of the i-th row and the j-th column and the 
zero otherwise. The indices i, j, k and I run at most from 1 to n. Only additional 
constraints on the indices are indicated. The subscript p runs from 1 to s, the 
subscript q runs from 1 to s'. The summation convention over repeated indices p 
and q is used unless otherwise stated. The number s is in the range 0, . . . , n — 1. 
In the case s = we assume 7 = 0, and all terms with the subscript p should be 
omitted from consideration. The value s' {s' < s) is defined below. 

Thus, the basis elements e.y ~ Efj, i < j, and fp r-^ YlilpiE^ satisfy the 
commutation relations [eij,eiiji\ = Si'j^ij' — ^ij'ei'j, [fp^^ij] = ilpi — 7pj)^ijj where 
6ij is the Kronecker delta. 

The Lie algebra t^(n) can be considered as the Lie algebra of the Lie subgroup 
T^{n) = {B e T{n) \ 3 £p e ¥ : bu = e^p^^p} of the Lie group T{n) of non-singular 
upper triangular nxn matrices. 

Below Al^'^l^, where ii ^ 12, ji ^ j2, denotes the submatrix of ^ 

matrix A = {aij). The conjugate value of k with respect to n is denoted by x, 
i.e., X = n — k + 1. The standard notation |^| = det A is used. 
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Proposition 2. Up to the equivalence relation on algebra parameters, the follow- 
ing conditions can be assumed satisfied for some s' € {0, . . . ,min(s, [n/2])} and 
kq, q = 1, . . . ,s' , 1 ^ /ci < A;2 < • • • < fc^/ ^ [n/2]: 

Iqk = lq>c, k < kq, -fq^^ - -fqj,^ = 1, 7pfc^ = Jp^^, p ^ q, q = I, . . . , s' , 

Ipk = Ipx, P> s', k = l,..., [n/2]. 

We will say that the parameter matrix 7 has a reduced form if it satisfies the 
conditions of Proposition [2l 

Theorem 1. Let the parameter matrix j have a reduced form. A basis o/Inv(ty(n)) 
is formed by the expressions 

s' 

[£]il[ n l^i^^r^N G {1, • • • , [n/2]} \ {fci, . . . , 

9=1 



-1)^+1 
fc=l \'^>t,n\ 



k<i<3 



pl,k r-l.A; 



p = s' + l,...,s, 

where x := n — k + 1, ii ^ «2; Ji ^ 32, denotes the matrix (eij)*C7'"''*|, and 

k 

aqk--= - ^(t?^- 7gfc')- 



fc'=i 



We use the short 'non-symmetrized' form for basis invariants, where it is uni- 
formly assumed that in all monomials elements of is placed before (or after) 
elements of f !^ n ■ 



6 Conclusion 

The main advantage of the proposed method is in that it is purely algebraic. Un- 
like the conventional infinitesimal method, it eliminates the need to solve systems 
of partial differential equations, replaced in our approach by the construction of 
the matrix B(Q) of inner automorphisms and by excluding the parameters Q from 
the fundamental lifted invariant I = x ■ B{6) in some way. 

The efficient exploitation of the method imposes certain constraints on the 
choice of bases of the Lie algebras. See, e.g.. Proposition [2] and Theorem [H That 
then automatically yields simpler expressions for the invariants. In some cases 
the simplification is considerable. 

Possibilities on the usage of the approach and directions for further investiga- 
tion were outlined in our previous papers [3-7]. Recently advantages of the moving 
frames approach for computation of generalized Casimir operators were demon- 
strated in [19] with a new series of solvable Lie algebras. The problem on optimal 
ways of applications of this approach to unsolvable Lie algebras is still open. 
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